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Abstract 

Double-lepton polarization asymmetries in the B — > £ + £~j decay are calculated 
using the most general, model independent form of the effective Hamiltonian including 
all possible forms of the interaction. The dependencies of the asymmetries on new 
Wilson coefficients are investigated. The detectability the averaged double-lepton 
polarization asymmetries at LHC is also discussed. 
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1 Introduction 



Rare radiative leptonic B s ^) — > £ + £~l decays are induced by the flavor-changing neutral 
current transitions b — > s(d). In the standard model (SM) such processes are described by 
the penguin and box diagrams and have branching ratios 10 -8 — 10~ 15 (see for example 
[1]). These rare decays can not be observed at the running machines such as Tevatron, 
BaBar and Belle, but the B s ^ — > ji + fj,~ and B s ^ — > /i + /i~7 decays are already studied 
[2]. Moreover, the transition form factors and many experimental observables such as, 
the branching ratio, photon energy, dilepton mass spectra and charge asymmetries are 
calculated for the B s u) —> £ + £~l decays in [3-9]. At the same time B s r,n — > £ + £~^i decays 
might be sensitive to the new physics beyond the SM. New physics effects in these decays can 
appear in two different ways: either through the new operators in the effective Hamiltonian 
which are absent in the SM, or through new contributions to the Wilson coefficients existing 
in the SM. One efficient way for precise determination of the SM parameters and looking for 
new physics beyond the SM is studying the lepton polarization effects. It has been pointed 
out in [10] that some of the single lepton polarization asymmetries might be too small to be 
observed and might not provide sufficient number of observables in checking the structure 
of the effective Hamiltonian. In need of more observables, in [10], the maximum number 
of independent observables have been constructed by considering the situation where both 
lepton polarizations are simultaneously measured. 

In the present work, we analyze the possibility of searching for new physics in the 
B — > £ + £~7 decay by studying the double-lepton polarization asymmetries, using the most 
general, model independent form of the effective Hamiltonian including all possible inter- 
actions. 

The work is organized as follows. In section 2, the matrix element for the B s — > £ + £~7 
is obtained, using the general, model independent form of the effective Hamiltonian. In 
section 3, we calculate the double-lepton polarization asymmetries. Section for is devoted 
to the numerical analysis, discussions and conclusions. 

2 Matrix element for the B q — > £ + £~j decay 

In this section we derive the matrix element for the B — > using the general, model 

independent form of the effective Hamiltonian. The matrix element for the process B — > 
£ + £~7 can be obtained from that of the purely leptonic B — > £ + £~ decay. At inclusive 
level the process B — > £ + £~ is described by b — > q£ + £~ transition. The effective b — * q£ + £~ 
transition can be written in terms of twelve model independent four-Fermi interactions in 
the following form [11]: 



= ^V tq V t l\c SL qi°J~ 2 Lbi-fi + C BR qia^Rb^H 




+ ClI QLi»b L l Ll n L + cf R q LliX b L l Rl n R + c RL q Rltl b R l Ll n L 
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where Cx are the coefficients of the four-Fermi interactions and 

2 ' 2 

The terms with coefficients Csl and Cbr which describe penguin contributions correspond 
to —2m s Cj^ and —Im^C^^ in the SM, respectively. The next four terms in this expression 
are the vector interactions. The interaction terms containing Cr°r and Cr°R in the SM have 
the form - C 10 and + C w , respectively. Inspired by this c£ knd C& will be 
written as 

^LR — ^9 + WO + <^LR , 

where Cll and Clr describe contributions from new physics. The terms with coefficients 
Clrlr, Crllr, Clrrl and Crlrl describe the scalar type interactions. The last two terms 
in Eq. (1) with the coefficients Ct and Cte describe the tensor type interactions. 

Having presented the general form of the effective Hamiltonian the next problem is the 
calculation of the matrix element of the B q — > £ + £^'~f decay. This matrix element can be 
written as the sum of the structure-dependent and inner Bremsstrahlung parts 

M = M S d + M ib . (2) 

The matrix element for the structure-dependent part Msd can be obtained by sandwiching 
the effective Hamiltonian between initial B and final photon states, i.e., (7 \H e ff\ B). It 
follows from Eq. (1) that in order to calculate A4sd, the following matrix elements are 
needed 



( 7 |S7„(1=F75)&|£) , 
{l\sa^ v b\B) , 

( 7 |s(l=F75)&|S) • 
The first two of the matrix elements in Eq. (3) are defined as [3, 12, 13] 



(3) 



( 7 (fc) 1^(1 =F 75)6| B(p B )) = ^{e, uXa e*»q x k°g(q 2 ) ± i [e*»(kq) - (e*q)k»]f(q 2 )} , (4) 



(l(k) |g<v&| B(p B )) 



ml 



Ge* x k a + He* x q u + N(e*q)q x k <T ] , 



(5) 



respectively, where e* and k are the four vector polarization and momentum of the photon, 
respectively, q is the momentum transfer, ps is the momentum of the B meson and g(q 2 ), 
f{q 2 ), G(q 2 ), H(q 2 ) and N(q 2 ) are the B — > 7 transition form factors. The matrix element 
( 7 (fc) |sc^i/ 7 5&| B{p B )) can be obtained from Eq. (5) using the identity 



_ 6 , 



2 



The matrix elements ("/(k) =f 7s)&| B(pb)) and (7 \sia tiv q u b\ B) can be calculated by 
contracting both sides of the Eqs. (4) and (5) with q^ and q u , respectively. We get then 

( 7 (fc)|s(l=F75)&|£(p B )>=0, (6) 
(7 \sia^q»b\ B) = -^i e^e^G . (7) 



m 



The matrix element (7 \sia^ v q u {l + 75)6! B) can be written in terms of the form factors that 
are calculated in framework of the QCD sum rules [3] and light front model [4] as follows 

(^\sia^(l+^)b\B) = A^{e^ a s a *qPk° gi {q 2 ) +i\e* (1 (qk) - (e*q)kjf 1 (q 2 )} . (8) 



It should be noted that these form factors were calculated in framework of the light-front 
model in [13]. Eqs. (5), (7) and (8) allow us to express G, H and N in terms of the form 
factors gi and f\. Using Eqs. (4)-(8), M.sd can be expressed as 

M SD = ^rV th V: q A^{h^l- l5 )t \A 1 e, ua ^q a k^ + tA 2 (e;(kq)-(e*q)k fl ) 

+ 1^(1 + 7sK [B^^q^ + i B 2 (e;(kq) - (e*q)k^) 
+ ie^pta^e \Ge* a k p + He^q 13 + N{e*q)q a k^ 



(9) 



+ ila 



Ci(£^r - eT) + H^cf - e*"<f) + N 1 (e*q){q fM k u - q v k») 



where 



A 2 
B 1 
B, 



j 2 (C BR + C 8L )gi + (C% + C RL )g , 

-^2~{CbR - CsLjfl + (C^l - C RI )jf , 

j 2 (C BR + C 8L )gi + (C& + C RR )g , 

~i {pBR — C S l) fl + (pLR — C RR j f , 



G = AC T gi , 
N=-±C T \{h + 9l ) 

H = N(qk) , 
G\ = —8Cte9i , 

N 1 = 8C TE \(f 1 + g 1 ) 
q A 

H x = N^qk) . 
For the inner Bremsstrahlung part we get 

M IB = -^-V tb V t * q ef B t^Fl 



(10) 
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2pik 2p 2 k 
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where we have used 



(0\sj^ 5 b\B) = - ifsPBn , 
(0|s<v(l + 75 )&|£) = 0, 

and conservation of the vector current. The functions F and Fi are defined as follows 
F = 2mA C L ° R — C L °[ + C RL — C RR ) H — (Clrlr — Crllr — Clrrl + Crlrl) , 

TIT? 1 

F\ = — — (Clrlr — Crllr + Clrrl — Crlrl] ■ (12) 



3 Double— lepton polarization asymmetries in B q — > 
£ + £~j decay 

In the present section we calculate the double-lepton polarization asymmetries, i.e., when 
polarizations of both leptons are taken into account. In order to calculate the double 
lepton polarization asymmetries we define the following orthogonal unit vectors sf in the 
rest frame of £ ± 



S L 



S L 



0,e L " 
0, 

0, et 

0,e + 
0,e£ 



P- 



\P-\J 
Pa x p- 



I ' 

( 0, . 

\ |PAXP-| 

(0, x eZ~) 
P+ \ 



0, 



\p+\) 

Pa x p+ 
\Pa x p+| 



= (0, x e/) 



(13) 



where p± and A; are the three-momenta of the leptons and photon in the center of mass 
frame (CM) of t~ i + system, respectively. Transformation of unit vectors from the rest 
frame of the leptons to CM frame of leptons can be accomplished by the Lorentz boost. 
Boosting of the longitudinal unit vectors yields 



( 



cm \mp 'me 



(14) 



where p + = — p_, Eg is the energy of the lepton in the CM frame and me is its mass. The 
unit vectors s^, are unchanged under Lorentz transformation. 

Having obtained necessary expressions, we now define the double-polarization asymme- 
tries as follows [10]: 
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dTisi ,sf) dr(-Si ,s" 



dq 2 



dq 2 



-( 



dT(s,r,-a+) dT(-s-,-s+) 



dq 2 



dq 2 



'drfa ,57) , dT(-Si ,s 



dq 2 



+ 



dq 2 



dT^i ,-57) , rfr(-s; ,-57)- 



(15) 



+ 



dq 2 



where, the first subindex i represents lepton and the second one antilepton. Using this 
definition of P^-, and after lengthy calculations, for the nine double-lepton polarization 
asymmetries we get 



Pll = 



—m B s 2 {l - s)(l - v 2 )lm[(A* 2 + B* 2 )HA 

run 



+ 16s(l - s)Re[v 2 G*H - (1 - 2v 2 )G* 1 H 1 ] 



+ 16s z 



v 2 \H\ 2 - (1 - 2v 2 ) \H X \ 



+ \fK{^ ~ s?{Ti +X 4 ) - [2a + (1 + s 2 )v 2 ]l 3 + [2a - (1 + s 2 )v 2 ]l,} \F\ 2 

+ \fW B { - (1 - s) 2 Zi + t; 2 [l + s(2 - s + 2sv 2 )}l 3 

- (1 - s)[(l - s)J 4 + 2su(l - w 2 )X 5 ] + w 2 [l - s(2 + s - 2sv 2 )]l 6 } \F x \ 2 

- 4f B m 2 B sv[(l - s)vls + (1 + s)l 9 ]Re[F*H] 

- Af B m%s[{\ - s)v 2 l 8 - (1 - s - 2v 2 )l 9 ]lm[F*H 1 ] 

+ /sw|s[8(l - s + 2sV) + m|(l - s)(4 + v 2 - 4s + 3sv 2 + 2sv 4 )J 8 

+ m|(l - s)(3 - 3s - 2v 2 + 4sv 2 )T 9 ]]lm[F*iV 1 ] 

+ f B [8(1 - s - 2sv 2 ) + m|(l - s)(4 + w 2 - 4s + 3sv 2 - 2sw 4 )J 8 

- m|(l - s)(l - s + 4v 2 - 2sV)J 9 ]]lm[F*G' 1 ] 

+ f B [ - 8(1 + s) + m|(l - s)(4 - v 2 - 4s + 3st; 2 )J 8 + m|(l - s)(l + s - 4-y 2 )J 9 ]] Re[F*G] 

+ / B m|s[ - 8(1 + s) + m|(l - s)(3 - s)t; 2 J 8 + m|(l - s 2 )J 9 ]Re[F*iV] 
1 



f B m B 4s(l + v 2 - s + 3sv 2 ) + m B s(l - s)(4 - 4s - 3v 2 + 7sv 2 + u 4 - sv 4 )l 8 



- m 2 B s(l - s)(l - s + v 2 + 7sv 2 - 4sv 4 )J 9 ]]Rep; + B^F^ 
1 



f B m B s 8(1 + s)v 2 + m 2 B (l - s)(2 - 2s - 2w 2 + 2sV + + sw 4 )J s 



2m 

m|(l - s> 2 T 9 ]]Re[(A* + B{)F] 
|(l-s) 2 (l-3, 2 )(|G| 2 + |G 1 | 2 ) 



+ ^-m B s(l - s) 2 (l 
3m£ 



2 ){m|slm[(^ + S^JVx] - 2Im[(A* + B^G,}} 



+ —m B s(l - s) 2 (l - t; 2 )Re[K + B{)G] 
6rri£ 



' mls^l-sfil-v^MAlB^ A* 2 B 2 ] 



3mj 

2 2 

rr 

3 



- -m B s{l -s) 2 (l + 3t; 2 )(|A 1 | 2 + \A 2 \ 2 + ^ + \B 2 \ 2 ) 



3 
4 



+ -m 2 B s(l - s) 2 (1 - 2v 2 )Re[G* 1 N 1 ] - v 2 Re[G*N] 



3 



- i<s 2 (l - s) 2 v 2 \N\ 2 - (1 - 2v 2 ) {mi 2 \ , (16) 



ln A 



^2f 2 m B m e VI(l - syimlF.F*}^ +1 4 ) 

+ f B m%Vi(l - s 2 )^ 2 Im[At(F 1 + F) + E*^ - F)]J 7 

+ 8f B m 2 B m e VI(l - s)v 2 lm[F*G]l 7 

+ 7im B V§{l - s) VRe[(A* -A* 2 + B{ + Bfid] 

- irm B V§(l - s)Vlm[(A* - A* 2 - B{ - B* 2 )G] 

+ 27rm B VP{l - s)w 2 [Re[(A* + B*)^] - Im[(A* - £*)#] 

- 47r/ B m B Vi(l - s)(l - v / T^ 2 )Im[(^ - S*)Fi + + £ 2 *)F] j , (17) 

P^vl = ^ 1 2/ 2 m B m^v / I(l - s)u 2 Im[FiF*](J 2 + J 4 ) 

+ /Bm|>/i(l - s 2 )^ 2 Im[At(F 1 - F) + ^ (F + F)]I 7 

+ 8f B m 2 B m e y/§(l ~ s)v 2 lm[F*G]l 7 

+ ttwibV^I - s) VRe[(A* + A* 2 + B{ - 

+ titob^I - s) Vlm[(A* + ^ - 5* + 5 2 *)G] 

+ 2*7^^(1 - s)v 2 [Re[(Al + + Im[(A* - £*)#]" 

+ 47r/ B m B v / I(l - s)(l - v / T^ 2 )Im[(A 2 ' - 5*)^ - (Aj; + £ 2 *)F] j , 
P LT = i| - -^f 2 B m B m e (l - s)v\(l - s) |F X | 2 + (1 + s) \F\ 2 }(1 2 +1 4 ) 



+ 8/ s m|m^(l - s^Imfi^ifi] + Re[F*H})l 7 

+ f B m 3 B VI(l - s) 2 vRe[(A* 2 + B^F^l, 

+ ^f B m B \Ts{l -s){\- VI - v 2 )Re[(A* 2 - B* 2 )F] 

- ^nf B m B VI[l - s(l - 2v 2 )](l - v / T^2)Re[(A* - 5*)^i] 

- A7rm,(l-s) 2 t;(|G' 1 | 2 + |G' 1 | 2 

- &irih e V§(l - s)uRe[GJifi + 



(18) 
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+ nm B V§(l - s) 2 vlm[(Al - A* 2 - B{ - Bfid] 

+ Trm B V§(l - s) 2 v[Re[(Al -A* 2 + B{ + B* 2 )G] + 2m B mA\e[A\A 2 - B{B 2 ] 
+ 2^5^(1 - s)v[lm[(Al - B^Hi] + Re[(A* + B{)H\ 









V 




32 


+ 


\/lv 




32 


+ 








Ptl = 


H 



irme(l - s)(l - \fl - ^Im^d] 
7rm £ (l - y/l -v 2 )Re[F*G}^ , 



+ 8f B m 2 B m e V§(l - s)w(lm[F*iJi] + Re[F*H])l 7 

+ f B m%V§(l - s) 2 vRe[(A* 2 + B* 2 )FAX 7 

- ^nf B m B VI(l - s)(l - Vl - v 2 )Re[(A* 2 - B* 2 )F] 

+ -nf B m B VI[l - s(l - 2t; 2 )](l - v / I^)Re[(^ - B*)^] 



- _7rm*(l-s) 2 7;(|Gi| 2 + |Gi| 

- 87rm*>/I(l - s)uRe[GJifi + 

- Ttm B \fi{\ - s) 2 vlm[(Al +A* 2 - B{ + Bfid] 

+ nm B V§(l - s) 2 v[Re[(Al + A* 2 + B{ - B* 2 )G] - 2m B m e Re\A* 1 A 2 - B{B 2 

- 27rm B V¥(l - s)v[lm[(Al - B{)H]\ - Re[(A\ + B{)h\ 

- -irf B m B V§{l + s)(l - v / T^ 2 )Re[K + B{)F] 

v 

+ -— 7rm<(l - s){\ - VT^jlmlF^d] 



32 1 

+ -y^7vm e (l - VT^)Re[F*G] \ , 



NT 



a' 



- 16s(l - s)uRe[G*i/i + G\H] 



- 16m B m e s(l - s)vlm[(A* 2 + B*)H] 



16, , 



2w - (1 -v 2 )ln 



1 + w 
1 -v 



Re[F*Gi 



+ 2$ B m\m t {\ - sfvImiAKF, - F) + 5*^ + F)](J 8 - J 9 ) 

- 32s 2 vlm[H*H] 

- f 2 B m B slm[F*F]{v[3 + v 2 - s(l - v 2 )]l 3 - (1 - v 2 )[(l - s)l b + (1 + s)J 6 ]} 

- 4/ s m B (l -s)[s- 2m 2 (l + s)]vIm[F*iV]J 8 



(19) 



(20) 
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+ 16/bs 



2v - (1 -v 2 )ln 



1 + v 



lm[F*G] 



- 2j B m\m i {l - s 2 )vlm[(A* 2 + 5 2 *)F](T 8 - J 9 
+ 4/ B m|suRe[F*ifi] [(1 - s)l 8 + (1 + s)Xg~ 
+ 4/ B m|svIm[F 1 *#] [(1 - s)v 2 I 8 + (1 - s + 2sv 2 )l 9 



16, 2 . 
+ -^fBm B s 



2v - In 



1 - v 



Re[F*JVi 



~fBm B m e (l - s) 



2v- (1 -v 2 )ln 



1 + v 
l-v 



Im[(A* - S*)F!] 



16 



m B m^(l - s) 2 uRep* + A* 2 + 5* - + m|s(A* - fl^Ni] 



- -m B (l - s) 2 vlm[2m e {Al + A* - B{ + B* 2 )G - m B s{ A iB\ + A* 2 B 2 )} 



s(l - s) 2 v(Re[G*7Vi + GJiV + m|siV*7V] + m B mJm[(A* + Bi 



• N = — | - 16s(l - s)wRe[G*// 1 + G*if] 
- 16m s m^(l - s)vIm[(A* + B*)H] 



16 , A 



2u - (1 - w 2 ) In 



1 + v 



Re[F*Gi] 

+ 2f B m B m l {\ - s) 2 uIm[AJ(Fi + F) + B*^ - F)](J 8 - J 9 ) 

- 32s 2 vIm[ J H 1 *iJ] 

- / 2 TO B sIm[F*F]{v[3 + v 2 - s(l - v 2 )]J 3 - (1 - v 2 )[(l - s)X 5 + (1 + s)X 6 ]} 

- 4/ s m B (l - s)[s - 2m 2 (l + s)]vlm[F* N]I 8 

l + v\l 

x Im[F*G] 



+ 16/bs 



2v - (1 - v 2 )ln 



- 2f B m B m e (l - s 2 )vlm[(A* 2 + 5 2 *)F](J 8 - J 9 ) 
+ AfsmlsvMF*^] [(1 - s) J 8 + (1 + s)l 9 

+ Af B m 2 B sv\m\FlH} [(1 - s)v 2 J 8 + (1 - s + 2sv 2 ) J 9 



16 



2v - In 



RefF*^ 



+ —f B m B m e (l - s) 

v 

16 



2v- (1 -v 2 )ln 



1 + v 
l-v 



lm[(A* 2 - B* 2 )F 1 ] 



- —m B m e (l - s) 2 uRep* -A* + B{ + B* 2 )G 1 - m B s(A\ - B\)N^ 

+ -m B (l - s) 2 vlm[2m £ {Al - At, - B{ - B* 2 )G - m B s{ A * B \ + A l B 2)\ 
3 

+ |m|s(l - s) 2 v(Re[G*iV 1 + G*iV + m|s7V*7V] + m B m e lm[(A* 2 + 5*)iV]) 
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Pnn = -^<- Wm B m e s(l - s)lm[(A* 2 + B^H,] 



+ lQsv 2 [(1 - s)Re[G*H] + s \H\ 2 
16s(l - s)Re[G* l H l - v 2 G*H] 



+ f 2 B m%s[(l + v 2 )l 3 -(l-v 2 )l 6 

- f 2 B m B sv{v[2 - s(l - v 2 )l 3 - (1 - v 2 )[(l - s)l 5 + svl 6 ]} \Fi\ 

- Af B m 2 B sv 2 Re[F*H] [(1 - s)T 8 + (1 + s)J { 
+ AfBmlsIm^Hi] [(1 - s)v 2 l 8 + (1 - s + 2sv 2 )l 9 



+ —fBm B m e s 

v 

16 , A 
+ — /bs 



2u - (1 - v 2 ) In 



1 + u 



2w - (1 -v 2 )ln 



1 + w 



+ ~fBm B s 
v 

- —fBm 2 B s 

V 

A 



l-v 

Av — (3 — s — v 2 + s-u 2 ) In 
1 + v 



ReUl + B^F,] 

(lm[F*Gi] - Re[F*G]) 
1 + v y 



Im^iVi 



2w - In 



l-v 



Re[F*iV] 



+ -m|s(l - s) V^RefA*^ + - G*N] - m|s |iV| 2 ) 
+ -m 2 B s(l - s) 2 (3 - 2v 2 )(2Re[GlN 1 ] + mislay 



(23) 



P TT = _|i6m B 7M(l - s)Im[(A* + 
- 16sw 2 [(l -s)Re[G*i/] + s|#| 2 " 
^/M{-(l-s) 2 (l-^i + [(l-*) 



3sv 2 (2 



2sV]J 3 



- (1 - v 2 ){l - s) 2 l 4 - 2s(l - s)v(l - v 2 )X 5 + (1 - v 2 )[l -s(2-s + 2sv 2 )]l 6 } |F X 

+ " (1 " S) 2 (l - t; 2 )Xi + [1 - v 2 - As + s 2 {\ - v 2 )]l 3 

- (1 - W 2 )(l - sfl A + (1 - t; 2 )(l - s 2 )T 6 } |F| 2 

- Af B m 3 B m e (l - s) 2 Re[(Al + B*)F](T 8 - J 9 ) 
+ Af B m 2 B sv 2 Re[F*H][(l - s)l 8 + (1 + s)l 9 ] 

- 4/ B m|sIm[F 1 *ffi][(l - s)t; 2 J 8 + (1 - s + 2st; 2 )J 9 ] 

- f B m B m e (l - s)Re[(^ + #2)^1] [§ - m|(4 - v 2 - 4s + 5sv 2 )J 8 
+ m|(3- 3s + 4sv 2 )J 9 

+ / B m|sRe[F*iV] [8(1 + s) - m|(l - s)(3 - s)w 2 J 8 - m|(l - s 2 )T 9 

- /smlsImfF*^!] [8(1 -s + 2sv 2 ) - m|(l -s)(A-v 2 -As + 5sv 2 - 2sv A )X 8 
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+ m|(l - s)(l - s - 2v 2 )l 9 
- / B Im[F*G' 1 ] [8(1 -s + 2sv 2 ) - m|(l - s)(4 - hv 2 - As + 9sv 2 - 2sv 4 )l 8 
+ m|(l - s) (3 - 3s - Av 2 + 6sv 2 )X 9 



+ f B Re[F*G}\8(l + 



m 



B 



{I - s){A - As - 3v 2 + sv 2 )T 8 



- m|(l-s)(5-3s-4t; 2 )J 9 ] 
+ ym B m,(l - s) 2 [2Im[(^ - B^G^ - 2Re[{A\ + 
+ m B m,(|A 1 | 2 + |A 2 | 2 + | J B 1 | 2 + | J B 1 



+ ^ 2 (1-5) 2 (|G 1 | 2 + |G| 2 ) 



3s 

+ -m|(l - s) 2 (sRe[A*5! + AJB 2 ] + m B m e lm[{A* 2 + B* 2 



+ -m 4 B s 2 (l - s) 2 [(l - 2v 2 ) m 2 + v 2 \Nf 



+ - 3 m B 



s(l - s) 2 Re[(l - 2v 2 )G\N l + v 2 G*N] j , 



where, 



(24) 



A = lQm B rhfXl - s) 2 (lm[(A; + B^d] - Re[(A* + B\)G - m B m t {A\B x + A\B X 
+ 48m B m^s(l - s)lm[(A* 2 + B*)^} 
- 8m B m e s(l - s) 2 Im[(A; + B*)^} 

+ \{l - s) 2 [4(3 - v 2 )( |G X | 2 + \G\ 2 ) + m%8(3 + v 2 )^ 2 + \A 2 \ 2 + \B,\ 2 + |B 2 | 2 ) 



+ 16sV 



(1 - s)Re[G*H}+ s\H\ 



+ 16s(3 - 2v 2 ) [(1 - s)Re[G^i] + s\Hi\ 
- -mls{l-s) 2 {Z-2v 2 )(me[G\N l }+m 2 B s\N 1 \ 2 ) 



-m, 



- s) V(2Re[G*iV] + m|s | TV | 2 ) 

- ^ |F| 2 {(1 - s)V(Jx + J 4 ) - (1 + s 2 + 2s^ 2 )T 3 - [1 - 5(4 - s - 2v 2 )]l 6 } 

+ \fW B \ p i I' { - (! - s)V(Zi + J 4 ) + [1 - s(2 - s - Av 2 + 2sv 2 - 2s^ 4 )]J 3 

- 2s(l - s>(l - v 2 )X 5 + [1 - s(2 - s + 2sw 2 - 2sv 4 )]X 6 ] 

- Af B m 2 B svRe[F*H][(l - s)vl 8 + (1 + s)J 9 ] 

- 4/ B m|sIm[F 1 *# 1 ][(l - s)w 2 J 8 + (3 - 2v 2 - 3s + 4sv 2 )J 9 ] 

+ 2f B m B m e Re[(Al + 5*)F] [8(1 + s) + m|(l - s> 2 J 8 + m|(l - s)(l - 3s)J 9 

- f B m B m e (l - s)Re[(A* + K[)Fi] f 8 + ~ 5s)v 2 J 8 + m|(3 - 3s + 4sv 2 )X 9 
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+ / B Im[F*Gi] [ - 24(1 -s + 2sv 2 ) + m|(l - s)(l + 3s - 6sw 2 )w 2 J 8 
+ to|(1 - s)(l - s - 2sv 2 )l 9 

+ / B Re[F*G] f - 24(1 + s) + m|(l - s)(l - 3s)v 2 X s - m 2 B (l - s)(l - 7s + 4sV)X 9 



+ / B m|sIm[W 



1 - s + 2sv 2 ) + m^(l - s)(3 + s - 2sv 2 )v 2 X f 



+ m|(l - s)(3 - 2v 2 - 3s + 4si; 2 )J 9 
+ / B TO|sRe[F*iV] f - 8(1 + s) + m|(l - s)(3 - s) V 2 J 8 + m|(l - s 2 )J 9 



In Eqs. (16)-(25), s = q 2 /m%, v — y 1 — 4to 2 /s is the lepton velocity with to^ = to^/to^, 
and 2j represent the following integrals 



where 



2"i = J i Ti{z)dz , 



^4 



^7 



(Pi • fc)(P2 • fc) 
-v2 



(Pi • &) 2 ' 



(P2 • ^) 2 ' 



, ft 



^8 



(Pi ' fc)(P2 • fc) 



(Pi ' ^) 2 ' 

-v2 



Pi • fc 



, ^3 



ft 



^9 



(Pi • fc)(p2 • k) ' 

1 



(Pi • k) 2 ' 
1 



Pi • k 



4 Numerical analysis and discussion 

We now proceed by presenting our numerical analysis for all possible double-lepton polar- 
izations. The values of the input parameters which have been used in the present work 



are: \V tb V t *\ 



0.0385, 



TO,, 



0.106 GeV, m T = 1.78 GeV, m b 



values of the Wilson coefficients we have used Cf M (m b 
r<SM 
°io 



-0.313, Cl M {m b ) 



4.8 GeV. For the SM 
4.344 and 

m b ) = —4.669. The magnitude of Cf M is quite well constrained from the b — > 57 
transition, and hence it is well established. Therefore the values of C B r and C S l are fixed 
by the relations Cbr = —2m b Cj^ and Csl = —2m s Cj^. It is well known that the Wilson 
coefficient C 9 M receives also long distance contributions which have their origin in the real 
cc intermediate states, i.e., J/ip, ip', ■ ■ ■ [14]. In the present work we restrict ourselves only 
to short distance contributions. 

In performing the numerical analysis, as is obvious from the expressions of given 
in Eqs. (16)-(24), we need to know the values of the new Wilson coefficients. During the 
numerical calculations, we will vary all new Wilson coefficients in the range 



°io 



< 



C x < 



°10 



and assume that all new Wilson coefficients are real. The experimental 
results on the branching ratio of the B — > K*(K)£ + £~ decays [15,16] and the bound on 
the branching ratio of B — > fi + fi~ [17] suggest that this is the right order of magnitude 
for the Wilson coefficients describing the vector and scalar interaction coefficients. But 
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present experimental results on the branching ratio of the B — > K*£ + £ and B — > K£ + £ 
decays impose stronger restrictions on some of the new Wilson coefficients. For example, 
-2 < C LL < 0, < C RL < 2.3, -1.5 <C T < 1.5 and -3.3 < C TE < 2.6, and all of the 



°io 



<C X < 



r<SM 
°io 



'16)-(24), explicit forms 



remaining Wilson coefficients vary in the region — 

In further numerical analysis, as can easily be seen from Eqs. 
of the form factors are needed, which are the main and most important parameters in the 
calculation of P^. These form factors are calculated in the framework of the QCD sum 
rules in [3, 12, 13] and their q 2 dependence can be represented, to a very good accuracy, in 
the following forms 

1 GeV „ ON 0.8 GeV 



9(<f) = ~, ; r~2 > /(<? 2 ) = 



o \ Z ' •> \ -± l / 9 \2' 

q \ q \ 



\ (5.6 GeV) 2 J \ (6.5 GeVf 

, „, 3.74 GeV 2 r . „, 0.68 GeV 2 

9i(<r) = -, - 9 , Mr) 



2 ' 



40.5 GeV 2 j \ 30 GeV 2 ) 

which we will be using in the numerical calculations. 

Numerical results are presented only for the B s — > £ + £~^ decay. It is clear that in 
the SU(3) limit the difference between the decay rates is attributed to the CKM matrix 
elements only, i.e., 

T(B d -> i + r 7 ) \v tb v t */ „ i 
r(B.->™- 7 ) |^| 2 20' 

It follows from the explicit expressions of the double-lepton polarization asymmetries 
that they depend on q 2 and the new Wilson coefficients. For this reason there may appear 
difficulties in studying the dependencies of the physical properties on both parameters at 
the same time. Hence, it is necessary to eliminate the dependence of Pij on one of these 
parameters. Here in the present work, we eliminate q 2 dependence of by performing 
integration over q 2 in the kinematically allowed region. The averaging of P^ over q 2 is 
defined as 

J4m ? 13 da 2 q 

(Pa) = J ( 26 ) 

The reason why we study the dependence of (Pij) on new Wilson coefficients is that in 
doing so we directly establish new physics beyond the SM, if the value of (Pij) turns out to 
be different compared to that predicted by the SM. 

In Figs. (l)-(8) we present the dependence of the averaged double-lepton polarization 
asymmetries on the new Wilson coefficients. From these figures we get the following results: 

• (Pll) exhibits very strong dependence on scalar, tensor interactions as well as on the 
vector interaction with coefficient Crl- When scalar interaction coefficients vary in 
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the region —4 < C sca \ aT < —0.8 the value of (Pll) is positive; in the region —0.8 < 
Cgcaiar < 0.8 it gets negative value and when it varies in the region 0.8 < C sca i ar < 4 
it again gets positive values. We should remind that in the SM {Pll) is negative and 
its magnitude is \{Pll)\ ~ 0.7. Similar situation holds for tensor interaction with 
the coefficient Cte-, as can easily be seen in Fig. (1). Our analysis shows that (Plt) 
depends more strongly on scalar and tensor interactions, as is the case for (Pll)- 
Therefore we can conclude that measurement of the sign and magnitude of (Pll) and 
(Plt) can give essential information about the existence of new physics beyond the 
SM. 

• (Ptl) is very sensitive to the existence of the scaler interaction only. More essential 
than that, (Ptl) changes its sign when scalar interaction coefficients vary in the 
allowed region. Such behavior can serve as a good test for establishing new type of 
scalar interaction. 

• (Ptt) ~ — (Pnn), and both are quite sensitive to the existence of tensor and scalar 
interactions, and to the vector interaction with coefficient Clr- In the presence of 
tensor, scalar and vector interactions the values of (Ptt) and (Pnn) can exceed the SM 
results 6, 5 and 2.5 times, respectively. Moreover, when scalar interaction coefficients 
Crllr(Crlrl) is negative, then the value of (Ptt) is positive (negative). When 
Crllr{Crlrl) is positive, (Ptt) becomes negative (positive). For the case when 
vector interaction coefficient Clr is negative, (Ptt) is also negative, while, when Clr 
changes sign and becomes positive, (Ptt) turns out to be positive as well. Therefore 
determination of the sign and magnitude of (Pnn) and (Ptt) can give unambiguous 
information about the existence of scalar and vector interactions. Departure of (Pnt) 
and (Ptn) from SM results is of not considerable importance, and hence we do not 
present their dependence on Cx- 

Depicted in Figs. (4)-(7) are the dependence of (Pij) on the new Wilson coefficients for 
the B s — > t + t~7 decay. Similar to the B s — > fi + fi~'~f decay, we observe that several of the 
double-lepton polarization asymmetries are very sensitive to the existence of new physics. 
More precisely, we can briefly summarize the results as follows: 

i) except Cll, (Pll) is very sensitive to the existence of all new Wilson coefficients. 

ii) (Plt) and (Ptl) exhibit strong dependence on all Wilson coefficients, except the 
vector interactions Cll and Crr. These quantities show very strong dependence, especially, 
on all types of scalar interactions, vector interactions with coefficients Clr and Crl, and 
tensor interaction with coefficient Ct- 

iii) (Ptt) = — (Pnn) are very sensitive to the existence of scalar interactions Crllr 
and Clrrl, when these coefficients are both positive. More important than that, (Ptt) = 
— (Pnn) change their sign when scalar coefficients vary in the region —3 > C sca i ar > +3. 

iv) (Pnt) ~ (Ptn) are both strongly dependent only on the tensor interaction with 
coefficient Cte- For all other new Wilson coefficients the values of (Pnt) ~ (Ptn) are very 
close to the SM prediction, i.e., to zero. Furthermore, when Cte is negative (positive), 
these quantities get positive (negative) values, with a considerable departure from the SM 
about 15%. Therefore, determination of the sign and and magnitude of (Ptn) and (Pnn) 
can give direct information solely about the existence of the tensor interaction. 
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v) Similar situation holds for the double-lepton polarization asymmetry (Pln) as well, 
i.e., (Pln) shows strong dependence only on the tensor interaction with coefficient Cte- 
When Cte is negative (positive), (Pln) gets negative (positive) values. Hence, departure 
from the SM prediction (in SM (Pln) ~ 0) can reach to 4%. Therefore analysis of (Pln) 
can serve as a good test for establishing the presence of tensor interactions. The values of 
(Pnl) is negligibly small (maximum departure from the SM is being about 1.5%) and for 
this reason we do not present its dependence on Cx- 

At the end of this section, let us discuss the possibility of measurement of the lepton 
polarization asymmetries in experiments. Experimentally, to measure an asymmetry (Pij) 
of the decay with the branching ratio B at no level, the required number of events (i.e., the 
number of BB pair) are determined by the following expression 

N 



Bs^Pij) 2 ' 

where si and s 2 are the efficiencies of the leptons. Efficiency of the /i lepton is practically 
equal to one, and typical values of the efficiency of the r lepton ranges from 50% to 90% 
for the various decay modes [18]. 

From the expression for M we see that, in order to obtain the double-lepton polarization 
asymmetries in B s — > £ + £~7 decays at 3a level, the minimum number of required events 
are (for the efficiency of r-lepton we take 0.5): 

• for the B s — > ^+^"7 decay 

~ 10 9 (Pll) , 

~7x 10 10 (P TT ) , 
~ 3 x 10 10 (P NN ) ~ (P LT ) ^ (Ptl) , 

which yields that, for detecting (Pln), (Pnl), (Pnt) and (Ptn), more than 10 13 BB 
pairs are required. 



M = < 



for B s — > t + t 7 decay 

~ 10 10 (Pll) , (Ptt) , (Pnn) , 



~ 3 x 10 11 (P LT ) = (Ptl) , 

> 10 13 (P LN ) , (Pnl) , (Ptn) , (Pnt) 



The number of BB pairs that will be produced at LHC is around ~ 10 12 . As a result 
of a comparison of this number of BB pairs with that of A/", we conclude that (Pll), 
(Ptt), (Pnn), (Ptl) and (P LT ) in B s -> ,u + /i~7 decay, and (P L l), (Ptt) and (P NN ) 
in B s — > r + r _ 7 decay can be detectable in future experiments at LHC. Note that for 
calculation of the branching ratio, we take its SM result, i.e., B(B S — > yU + /i~7) ~ 1.3 x 10 8 
and B(B S — > r + r _ 7) ~ 6 x 10 9 . In obtaining these values, minimal value for the photon 
energy is taken to be 50 MeV. 

In conclusion, we calculate nine double-lepton polarization asymmetries using the most 
general, model independent form of the effective Hamiltonian including all possible form 
of interactions. The sensitivity of the averaged double-lepton polarization asymmetries to 
the new Wilson coefficients are studied. Finally we discuss the possibility of experimental 
measurement of these double-lepton polarization asymmetries at LHC. 
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Figure captions 

Fig. (1) The dependence of the averaged double-lepton polarization asymmetry (Pll) on 
the new Wilson coefficients for the B s — > /i + /i"7 decay. 

Fig. (2) The same as in Fig. (1), but for the averaged double-lepton polarization asym- 
metry (Ptl)- 

Fig. (3) The same as in Fig. (1), but for the averaged double-lepton polarization asym- 
metry (Ptt)- 

Fig. (4) The same as in Fig. (1), but for the B s — > t + t - 7 decay. 

Fig. (5) The same as in Fig. (4), but for the averaged double-lepton polarization asym- 
metry (Plt)- 

Fig. (6) The same as in Fig. (4), but for the averaged double-lepton polarization asym- 
metry (P T l)- 

Fig. (7) The same as in Fig. (4), but for the averaged double-lepton polarization asym- 
metry (Ptt)- 
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